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Nuclear effective interactions are useful tools in astrophysical applications especially if one can 
guide the extrapolations to the extremes regions of isospin and density that are required to simulate 
dense, neutron-rich systems. Isospin extrapolations may be constrained in the laboratory by mea- 
suring the neutron skin thickness of a heavy nucleus, such as 208 Pb. Similarly, future observations 
of massive neutron stars will constrain the extrapolations to the high-density domain. In this con- 
tribution we introduce a new relativistic effective interaction that is simultaneously constrained by 
the properties of finite nuclei, their collective excitations, and neutron-star properties. By adjusting 
two of the empirical parameters of the theory, one can efficiently tune the neutron skin thickness 
of 208 Pb and the maximum neutron star mass. We illustrate this procedure in response to the 
recent interpretation of X-ray observations by Steiner, Lattimer, and Brown that suggests that the 
FSUGold effective interaction predicts neutron star radii that are too large and a maximum stellar 
mass that is too small. The new effective interaction is fitted to a neutron skin thickness in 208 Pb 
of only R n — R p = 0.16 fm and a moderately large maximum neutron star mass of 1.94 Mq. 

PACS numbers: 26.60.Kp,21.65.Cd,21.60.Jz 



I. INTRODUCTION 



Nuclear effective interactions provide a compact and efficient description of both the structure of finite nuclei and 
the equation of state — pressure as a function of density and temperature — of nuclear matter. Effective interactions 
are often fitted to well determined nuclear observables, such as binding energies and charge radii of closed shell nuclei. 
However, limiting the calibration procedure to these few properties leaves important parts of the effective interaction 
largely undetermined, particularly the isovector interaction. Failure to constrain the isovector interaction affects the 
properties of systems with large isospin asymmetry. In the case of finite nuclei, perhaps the best example is the 
neutron radius (or neutron skin) of heavy nuclei. Indeed, effective interactions calibrated to the binding energy and 
charge radii of closed-shell nuclei often predict significantly different neutron radii [HH]. In the case of infinite matter, 
it is the density dependence of the symmetry energy that remains largely unconstrained. These poorly constrained 
properties are very important in astrophysical systems with large neutron asymmetries, such as neutron stars [SHE]. 

In this article we develop an effective interaction suitable for the calculation of the ground-state properties of finite 
nuclei, their collective response, and the structure of neutron stars. Astrophysical predictions are notoriously difficult 
as they involve two large extrapolations away from the familiar landscape of stable nuclei. Whereas one involves 
isospin extrapolations from the domain of stable nuclei to the very neutron-rich systems present in neutron stars, the 
other one involves an extrapolation from the density of normal nuclei to the very high (and very low) densities that 
are encountered in neutron stars. To reduce the uncertainties associated with these extrapolations we impose two 
additional constraints — one experimental and one observational — on the effective interaction. Specifically, we extend 
the standard protocol of fitting to the binding energies and charge radii of finite nuclei by demanding that the new 
effective interaction also reproduces: (a) the neutron radius of 208 Pb and (b) the maximum mass of a neutron star. We 
have selected these two quantities over other possible choices because they can be determined in a model independent 
manner from experiment and observation. In particular, fixing the neutron radius in 208 Pb reduces significantly the 
uncertainties associated with extrapolations in isospin. The Lead Radius experiment (PREx) — which was successfully 
commissioned in March of 2010 — aims to determine the neutron radius of 208 Pb in a model independent fashion by 
using parity violating electron scattering [3 [TU] . To reduce the uncertainty associated with extrapolations to high 
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densities, we have selected to fix the maximum neutron star mass. Of course, the maximum neutron star mass selected 
for the fit must at least be equal to the largest — well measured — mass determined from observation [llj . 

In principle, other choices are possible to guide these extrapolations. For example, laboratory experiments with 
heavy ions have played a critical role in probing the nuclear equation of state. By tuning the energy of the colliding 
beams and the neutron-proton asymmetry, heavy-ion collisions probe vast regions of the phase diagram. In particular, 
isospin diffusion data may be used to constrain the density dependence of the symmetry energy [T2HT4] . This could be 
used instead of the neutron radius in 208 Pb. Moreover, experiments with very energetic heavy ions have compressed 
nuclear matter to densities in excess of four times nuclear matter saturation density |15) . The determination of the 
pressure of symmetric nuclear matter at these high densities may be used instead of the maximum neutron star mass. 
However, not only is the interpretation of heavy-ion data model dependent, but heavy-ion systems are much less 
neutron rich and often much hotter than neutron stars. This introduces large and uncontrolled uncertainties. 

Neutron star radii as a function of mass — the so-called Mass-vs-Radius relationship — provides critical information 
on the equation of state of neutron-rich matter at high densities. Yet, the simultaneous determination of both the 
mass and radius of a neutron star is difficult. However, tremendous advances in X-ray astronomy have produced large 
amounts of high quality data that have allowed the simultaneous determination of both masses, radii, and the underly- 
ing equation of state |161 117) . Unfortunately, neutron star radii inferred from X-ray observations of the luminosity and 
temperature are sensitive to models of the stellar atmosphere (see Ref. [18] and references therein). Moreover, mass 
and radius determination from X-ray burst (as in Ref. |16j ) may be hindered by systematic uncertainties [17] . For 
example Ozel, Baym, and Giiver inferred very small stellar radii from their recent analysis of three X-ray bursts 16J, 
while Steiner, Lattimer, and Brown extracted neutron star radii significantly larger (of the order of 12 km) [T7]. Still, 
we believe that such studies will eventually become instrumental in constraining the dense matter equation of state. 

In the past, we have studied correlation between the neutron skin of 208 Pb and the non-uniform solid crust of a 
neutron star [3J. For models with a stiff equation of state, namely, one where the pressure increases rapidly with 
density, it is energetically unfavorable to separate uniform nuclear matter into regions of high and low densities. 
Thus, models with a stiff equation of state predict low transition densities from non-uniform to uniform neutron-rich 
matter. The neutron skin thickness in 208 Pb also depends on the equation of state of neutron-rich matter: the stiffcr 
the equation of state the thicker the neutron skin. Thus, an inverse relationship was established: the thicker the 
neutron skin of a heavy nucleus, the lower the transition from non-uniform to uniform neutron-rich matter. This 
represents one of the many examples that show the utility of constraining the neutron skin thickness in 208 Pb — and 
consequently the effective interaction — for astrophysical applications. 

Many relativistic effective interactions exist. For example, the NL3 interaction [TS], [2U] provides an excellent 
description of the binding energy and charge radii of many nuclei throughout the periodic table. When extrapolated 
to the neutron-star domain, the NL3 effective interaction predicts a very stiff equation of state that generates both 
large stellar radii and a large maximum neutron star mass. On the other hand, the FSUGold (or FSU for short) 
interaction [3T] also provides a good description of closed shell nuclei, but predicts a significantly softer equation 
of state. Indeed, whereas NL3 predicts a maximum neutron star mass of 2.78 M , the FSU interaction predicts 
a limiting mass almost one full solar mass smaller (of 1.72 Mq). The softer equation of state emerges as one 
incorporates constraints from giant resonances and heavy-ion collisions. At present, we are not aware of any nuclear 
effective interaction that also incorporates astrophysical constraints in the calibration procedure. It is the aim of the 
present manuscript to obtain a new effective interaction that improves on both NL3 and FSU by incorporating some 
of the recent constraints suggested by Steiner, Lattimer, and Brown [17] . In particular, we will show that the EOS 
predicted by the new model generates a limiting neutron star mass intermediate between NL3 and FSU but stellar 
radii smaller than both. 

We have organized our article as follows. In Section [H] we describe details of the fitting procedure and present 
results for several ground-state properties of closed shell nuclei. Next, we show in Section [III] how the new effective 
interaction reproduces the properties of various nuclear collective modes. Predictions for the equation of state and 
for various properties of neutron stars are presented in Section |IV| Finally, we conclude in Section [V] 



II. RELATIVISTIC EFFECTIVE FIELD THEORY 



Our starting point will be the relativistic effective-field theory of Ref. [22] supplemented with an isoscalar-isovector 
coupling as introduced in Ref. [3] . The interacting Lagrangian density for this model is given by [3] 122] 
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The model contains an isodoublet nucleon field (4>) interacting via the exchange of two isoscalar mesons, the scalar 
sigma (0) and the vector omega (V M ), one isovector meson, the rho (b^ 1 ), and the photon (A 11 ). In addition to 
meson-nucleon interactions, the Lagrangian density includes scalar and vector self-interactions. (Note that while the 
original model allows for p-meson self-interactions [22j . their phenomenological impact has been documented to be 
small so they will not be considered in this contribution). The scalar self- interaction is responsible for reducing the 
compression modulus of nuclear matter from the unrealistically large value of K = 545 MeV [231 EI] all the way down 
to about K — 230 MeV. This latter value appears to be consistent with measurements of the isoscalar giant monopole 
resonance (GMR) in 208 Pb [25H2EI- 

Omega-meson self-interactions, as described by the parameter (, soften the equation of state at high density and 
can be tuned to reproduce the maximum mass of a neutron star. Indeed, Miiller and Serot found it was possible to 
build models with different values of £ that reproduce the same observed properties at normal nuclear densities, yet 
produced maximum neutron star masses that differ by almost one solar mass |22) . In particular, models with £ = 
predict large limiting masses of about 2.8M — even for models with a soft symmetry energy. In contrast, the nonlinear 
coupling constant A v was included to modify the density-dependence of the symmetry energy. Tuning A v provides 
a simple and efficient method of softening the symmetry energy without compromising the success of the model in 
reproducing well determined ground-state observables |21j . Here we will adjust A v by assuming a relatively small 
value for the neutron skin thickness of 208 Pb. Many properties of finite nuclei, such as the binding energies and charge 
radii of closed shell nuclei, are insensitive to the values of C and A v . Indeed, the NL3 effective interaction reproduces 
these observables without ever introducing these two empirical parameters. Our approach here is to accommodate 
newly proposed astrophysical constraints that suggest that, relative to the FSU interaction, the equation of state of 
stellar matter must be slightly softer at intermediate densities and stiffer at high densities [PT] . 

The new effective interaction is generated by adopting the following procedure. We start from the FSU parameter 
set [ST] as listed in Table [TJ This parameter set reproduces ground-state properties and collective excitations of closed 
shell nuclei. However, this interaction has a relatively large value of ( = 0.06. This reduces the pressure at high 
densities and generates a relatively small neutron star maximum mass of 1.72 M Q . Note that, to date, the only 
terrestrial constraint on the high-density component of the EOS comes from energetic heavy- ion collisions 15J. The 
FSU equation of state fits these data comfortable within the errors. Yet, the reported errors are large enough to 
accommodate slightly stiffer equations of state. As we aim for an interaction with a somewhat larger maximum mass 
of about 2.0 Mq, we reduce the value of ( from 0.06 to 0.03. Next, we refit the isoscalar parameters (g s , g v , k, and A) 
to maintain the saturation properties of symmetric nuclear matter at their FSU values, namely, a saturation density 
of 0.148 fm~ 3 (or a Fermi momentum of &f = 1-30 fm" 1 ), an energy per nucleon of E/A — —16.30 MeV, and an 
incompressibility coefficient of K = 230 MeV. 
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9t 


gi 


ai 


re 
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C 


A v 


NL3 
FSU 
IU-FSU 


508.194 
491.500 
491.500 


782.501 
782.500 
782.500 


763.000 
763.000 
763.000 


104.3871 
112.1996 
99.4266 


165.5854 
204.5469 
169.8349 


79.6000 
138.4701 
184.6877 


3.8599 
1.4203 
3.3808 


-0.015905 
+0.023762 
+0.000296 


0.00 
0.06 
0.03 


0.000 
0.030 
0.046 



TABLE I: Parameter sets for the three models discussed in the text. The parameter re and the meson masses m B , m v , and m p 
are all given in MeV. The nucleon mass has been fixed at M = 939 MeV in all the models. 



We then increase the isoscalar-isovector coupling constant A v from its FSU value of 0.03 to 0.046. This change 
softens the density dependence of the symmetry energy and reduces the neutron radius in 208 Pb. We aim for a modest 
value of the neutron skin thickness in 208 Pb of about R n — R p = 0.16 fm. This value, although smaller as compared 
to other relativistic mean field models, is close to those predicted by many non-relativistic models. Note that when 
we change A v , we also change g p in order to maintain the symmetry energy fixed at ~26 MeV at a Fermi momentum 
of fcp = 1.15 fm -1 . This procedure and the justification for it are described in detail in Ref. [3J. Increasing A v in this 
fashion reduces the density dependence of the symmetry energy. This generates a symmetry energy at low densities 



Model 


po (mi"*) 


so (MeV) 


K (MeV) 


J (MeV) 


L (MeV) 


NL3 
FSU 
IU-FSU 


0.148 
0.148 
0.155 


-16.24 
-16.30 
-16.40 


271.5 
230.0 
231.2 


37.29 
32.59 
31.30 


118.2 
60.5 
47.2 



TABLE II: Bulk parameters characterizing the behavior of infinite nuclear matter at saturation density p . The quantities e 
and Ko represent the binding energy per nucleon and incompressibility coefficient of symmetric nuclear matter, whereas J and 
L represent the energy and slope of the symmetry energy at saturation density. 
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that is larger than the one for FSU. As a result, protons near the surface of 208 Pb are pulled closer to the neutrons 
in an effort to minimize the neutron-proton asymmetry. This leads to an increase in the charge radius of 208 Pb. In 
order to maintain the charge radius in 208 Pb at its experimental value we must slightly increase the Fermi momentum 
of symmetric nuclear matter from 1.30 to 1.318 fm -1 . We then refit the parameters of the model to keep the other 
saturation properties of nuclear matter intact. Finally the scalar coupling g s is slightly tuned to improve the fit to 
the binding energies of closed shell nuclei. Note that we have not changed the scalar mass m s from its FSU value. 

The resulting parameter set — henceforth referred to as the Indiana University-Florida State University (IU-FSU) 
interaction — is listed in Table |TJ The resulting bulk properties of infinite nuclear matter have been collected in Table |H| 
In addition, predictions for several ground-state properties of closed shell nuclei are listed in Table III and compared 
against other theoretical models and experiment (when available). Finally, predictions for the charge and neutron 
densities of 208 Pb are displayed in Fig. [TJ Whereas significant differences can be observed in the prediction of the 
neutron densities, the difference among the models is small for the charge density. Note that although the parameters 
of the model will eventually be determined from an accurate calibration procedure, the new model reproduces rather 
well the charge radii and binding energies of closed shell nuclei. However, the IU-FSU interaction predicts the rather 
modest value of R n — R„ = 0.16 fm for the, as yet unknown, neutron skin thickness in 208 Pb. 



Nucleus 


Observable 


Experiment 


NL3 


FSU 


IU-FSU 


40 Ca 


B/A (MeV) 


8.55 


8.54 


8.54 


8.53 




Rch (fm) 


3.45 


3.46 


3.42 


3.41 




R n -R p (fm) 




-0.05 


-0.05 


-0.05 


48 Ca 


B/A (MeV) 


8.67 


8.64 


8.58 


8.55 




Rch (fm) 


3.45 


3.46 


3.45 


3.44 




Rn-Rp (fm) 




0.23 


0.20 


0.17 


9u Zr 


B/A (MeV) 


8.71 


8.69 


8.68 


8.67 




Rch (fm) 


4.26 


4.26 


4.25 


4.23 




Rn-Rp (fm) 




0.11 


0.09 


0.07 


laa Sn 


B/A (MeV) 


8.36 


8.37 


8.34 


8.33 




Rch (fm) 




4.70 


4.71 


4.68 




Rn-Rp (fm) 




0.35 


0.27 


0.22 


2U8p b 


B/A (MeV) 


7.87 


7.88 


7.89 


7.89 




Rch (fm) 


5.50 


5.51 


5.52 


5.48 




R n -R p (fm) 




0.28 


0.21 


0.16 



TABLE III: Experimental data for the binding energy per nucleon and charge radii for several doubly magic nuclei. Results are 
presented for the three models employed in the text alongside their predictions for the neutron skin thickness of these nuclei. 



III. GIANT RESONANCES 



At the time that the FSU effective interaction was calibrated, the only evidence in favor of a symmetry energy softer 
than NL3 came from nuclear collective excitations, particularly the mass dependence of the isoscalar giant monopole 
resonance (GMR). Whereas NL3 — with an incompressibility coefficient of K= 271 MeV was able to reproduce the 
centroid energy of the GMR in 208 Pb, it overestimated the GMR in 90 Zr. It was later confirmed that the success of the 
NL3 interaction in reproducing the breathing mode in 208 Pb was accidental, as it resulted from a combination of both 
a stiff EOS for symmetric nuclear matter and a stiff symmetry energy [3Dj. Indeed, with a relatively large neutron- 
proton asymmetry of ct= (N — Z)/A = 0.212, the GMR in 208 Pb probes the incompressibility of neutron-rich matter, 
rather than that of symmetric matter. Given that the incompressibility of neutron-rich matter softens with a [31] , 
NL3 could reproduce the GMR in 208 Pb by cancelling its stiff incompressibility coefficient K with a correspondingly 
stiff symmetry energy. This cancellation, however, is incomplete in 90 Zr because its nucleon asymmetry is almost 
twice as small as that of 208 Pb. Thus, to develop the FSU effective interaction we used the GMR in 90 Zr, rather than 
in 208 Pb, to fix the incompressibility coefficient of symmetric nuclear matter. Having fixed K, we used the GMR — 
and the isovector giant dipole resonance (IVGDR) — in 208 Pb to constrain the density dependence of the symmetry 
energy [3T]. As we shall see below (see Table IV), the newly proposed IU-FSU effective interaction continues to 
reproduce the centroid energy of these three modes. 

The distribution of both isoscalar monopole and isovector dipole strength were computed in a relativistic random- 
phase approximation (RPA). The first step in calculating the RPA response is the computation of the mean-field 
ground state in a self-consistent mean-field approximation. Once self-consistency is achieved, three important pieces 
of information become available: (i) the single-particle energies of the occupied orbitals, (ii) their single-particle wave 
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FIG. 1: (Color online) Model predictions for the charge and neutron densities of 208 Pb. The experimental charge density is 
from Ref. [29]. 



functions, and (iii) the self-consistent mean-field potential. This mean-field potential — without any modification — 
must then be used to generate the nucleon propagator; only then can one ensure the conservation of the vector current . 
The nucleon propagator, which is an essential building block of the uncorrelated polarization tensor, is computed non- 
spectrally to avoid any dependence on the artificial cutoffs and truncations that plague most spectral approaches. 
Moreover, a great merit of the non-spectral approach is that the continuum is treated exactly. Once the uncorrelated 
polarization insertion is computed, the correlated RPA polarization is obtained by iterating it to all orders with the 
consistent residual interaction. The distribution of strength is then obtained directly from the imaginary part of the 
polarization tensor. A detailed description of the relativistic RPA approach may be found in Ref. [MUM]- Results 
for the centroid energies of the giant monopole and dipole resonances are presented in Table |TV| for the three effective 
interactions considered in the text. Unlike the NL3 model, the results obtained with both the FSU and the IU-FSU 
effective interactions are consistent with experiment. 



Nucleus 


Observable 


Experiment 


NL3 


FSU 


IU-FSU 


2U8p b 

90 Zr 

208p b 


OMR (MeV) 
OMR (MeV) 
IVGDR (MeV) 


14.17 ±0.28 
17.89 ±0.20 
13.30 ±0.10 


14.32 
18.62 
12.70 


14.04 
17.98 
13.07 


14.17 
17.87 
13.24 



TABLE IV: Centroid energies for the GMR in 
data are extracted from Refs. [25] and 1331. 



3 Pb and 90 Zr, and the peak energy for the IVGDR in 208 Pb. Experimental 



IV. EQUATION OF STATE AND NEUTRON STAR STRUCTURE 

Although the full complexity of the quantum system can not be tackled exactly, the energy density of asymmetric 
nuclear matter may be computed in a mean-field (MF) approximation. In the MF approximation all the meson fields 
are replaced by their classical expectation values and their solution can be readily obtained by solving the classical 
Euler-Lagrange equations of motion. The sole remnant of quantum behavior is in the treatment of the nucleon field 
which emerges from a solution to the Dirac equation in the presence of appropriate scalar and vector potentials |23U34| . 
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Following standard mean-field practices, the energy density of the system is given by the following expression: 

gfaa) = \ k 2 Vk 2 + M* 2 dk+^ k 2 Vk 2 + M* 2 dk+\ (^) $g + ^ + 

Jo * Jo 2 V 9t J 6 24 

+ \ (S) W " + C 8 W " + * (if) Bl + 3AvW ° B2 ° ' (2) 

where p is the baryon density of the system and a — {p n ~Pp)/p is the neutron-proton asymmetry. Alternatively, the 
energy density may be written in terms of the proton and neutron Fermi momenta fcp and fcp. Note that the following 
definitions have been introduced: $o = 5s<A)j Wq= g v Vo, and BQ = g p bo- Also note that since the MF approximation 
is thermodynamically consistent, the pressure of the system (at zero temperature) may be obtained either from the 
energy-momentum tensor or from the energy density and its first derivative. That is, 

P{p,a)=p—-S. (3) 

The equation of state (pressure as a function of baryon density) for symmetric nuclear matter (fcp = fcp and B$ = 0) is 
displayed in Fig [2} By design, the FSU and IU-FSU interactions have (almost) the same incompressibility coefficient 
so they predict similar pressures at low to intermediate densities. However, at higher densities the IU-FSU interaction 
is stiffer because it has a smaller value of £ as compared to the FSU interaction. To date, the only terrestrial constraint 
on the high-density component of the EOS comes from energetic heavy-ion collisions. Both of these models predict 
equations of state that are consistent with the phenomenological flow analysis by Danielewicz, Lacey, and Lynch [15] . 
In contrast, the NL3 interaction, although enormously successful in reproducing ground-state properties of finite 
nuclei, predicts an equation of state that is significantly stiffer than the phenomenological extraction. 




FIG. 2: (Color online) The equation of state — pressure P vs baryon density — of symmetric nuclear matter. Here po is the 
density of nuclear matter at saturation and the shaded area represents the EOS extracted from the analysis of Ref. |15j . 



The equation of state of pure neutron matter (PNM) provides a powerful theoretical constraint on models of the 
effective interaction. By building on the universal behavior of dilute Fermi gases with an infinite scattering length, 
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significant progress has been made in constraining the equation of state of pure neutron matter [35] • To date, a 
variety of models using different neutron-neutron interactions and a variety of many-body techniques have been 
employed to compute the EOS of dilute neutron matter (see Ref. and references therein). In Fig. [3] we display the 
energy per particle of pure neutron matter for a variety of microscopic approaches |36| and for the three rclativistic 
effective interactions discussed in the text. In most relativistic descriptions the isovector interaction is modelled by the 
exchange of a single vector- isovector ( "rho" ) meson with its (Yukawa) coupling to the nucleon tuned to reproduce the 
symmetry energy at saturation density. Invariably, such a simple prescription generates a stiff symmetry energy. In 
turn, this yields an EOS for dilute neutron matter that is inconsistent with the model-independent results of Schwenk 
and Pethick (denoted by the shaded area in Fig. [3]) . The FSU isovector interaction improves on NL3 by adding an 
isoscalar- isovector mixing term (A v ) that softens the symmetry ener gy. Recall that such a softening is required to 
reproduce the centroid energies of various collective modes (see Table IV). Without any further adjustment, the FSU 
interaction is also consistent with the theoretical constraints (see solid blue line). Finally, the IU-FSU interaction — 
with an additional softening relative to FSU — remains within the theoretical "error bars" . 
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FIG. 3: (Color online) Energy per nucleon as a function of the Fermi momentum for pure neutron matter. 



We conclude this section by computing the mass-vs-radius (M-R) relationship of neutron stars. Given that the 
structure of neutron stars is sensitive to the EOS over an enormous range of densities, we must resort to various 
descriptions. For the non-uniform outer crust we employ the equation of state of Baym, Pethick, and Sutherland |37j . 
At densities of about a third of nuclear-matter saturation density, uniformity in the system is restored. For this 
(liquid-core) region we use an EOS generated from the relativistic models discussed in the text. Although significant 
progress has been made in simulating the exotic structure of the inner crust (i.e., the region between the outer crust 
and the liquid core) a detailed equation of state is still missing in this region. Hence, we resort to a fairly accurate 
polytropic EOS to interpolate between the solid crust and the uniform liquid core [351 |3j5] • The transition from the 
uniform liquid core to the solid crust is sensitive to the density dependence of the symmetry energy. The crust-to-core 
transition density p t is obtained from an RPA stability analysis to determine the onset of the instability to small 
amplitude density oscillations [3S]. With such an equation of state in hand, we can now proceed to integrate the 
Tolman-Oppenheimer-Volkoff equations to determine the structure of neutron stars. 

In Fig. [4] we display predictions for the equation of state of stellar matter (i.e., neutron-rich matter in beta equi- 
librium) for the three models discussed in the text, alongside the constraint extracted from Ref. [17] . For reference, 
also shown is the stiffest (P — 8) possible equation of state consistent with causality. As suggested in Ref. [T7], the 



FSU interaction appears slightly stiff at intermediate densities and too soft at high densities. The IU-FSU interaction 
successfully corrects both of these problems. Note, however, that the observational EOS can accommodate even stiffer 
equations of state. We stress that, if required, this can be implemented rather easily by reducing even further the 
value of (. 




500 1000 1500 2000 
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FIG. 4: (Color online) Equation of state — pressure vs energy density — of neutron-rich matter in beta equilibrium. The shaded 
region displays the observational constraint extracted from Ref. |17] . The solid black line (P = £) denotes the stiffest possible 
equation of state consistent with causality. 



In Fig. [5] we display predictions for the M-R relation and compare them against observational constraints extracted 
from the analyses of Refs. [TB] and [17] • The very stiff behavior of the EOS predicted by NL3 is immediately ruled 
out by both observational constraints. In the case of the FSU effective interaction, we recently identified a conflict 
when compared with the analysis by Ozel and collaborators [40] . Whereas the FSU model seems to generate an 
EOS that is consistent with observation, the predicted stellar radii were simply too large. Relative to the analysis 
by Steiner Lattimer, and Brown [17] . the FSU predictions overestimate the stellar radius — although not as much as 
suggested by [TB] — but underestimates the maximum neutron star mass. The IU-FSU interaction, with a softer EOS 
at intermediate densities and a stiffer one at high densities, is motivated in response to these findings. Indeed, the 
IU-FSU interaction predicts a maximum stellar mass of 1.94 M and a stellar radius of R — 12.49 km for a 1.4 Mq 
neutron star (see Table [v]) . These predictions are consistent with the 2a values extracted from X-ray observations by 
Steiner, Lattimer, and Brown [T7J. However, they are well outside the limits extracted by Ozel and collaborators 16 . 
Given that such extractions depend critically on the models used to simulate X-ray bursts, much work remains to 
be done to reconcile these two analyses. We close this section by listing in Table [V] predictions for several important 
neutron-star properties. In addition to the properties already discussed (such as masses and radii), the table includes 
the minimum density required for the onset of the direct Urea process and the minimum stellar mass that may cool 
down by the direct Urea process. Note that small neutron star radii and enhanced cooling are generally regarded as 
good indicators of a possible phase transition in the stellar core. 
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FIG. 5: (Color online) Mass-vs-Radius relation predicted by the three relativistic mean-field models discussed in the text. The 
observational data that suggest very small stellar radii represent la confidence contours for the three neutron stars reported 
in Ref. [16| . The two shaded areas that suggest larger radii are la and 2a contours extracted from the analysis of Ref. |17) . 



Observable 


NL3 


FSU 


IU-FSU 


Pt (fm" J ) 


0.052 


0.076 


0.087 


R (km) 


15.05 


12.66 


12.49 


M max (M ) 


2.78 


1.72 


1.94 


P U ^ (to" 3 ) 


0.21 


0.47 


0.61 


M Urca (M ) 


0.84 


1.30 


1.77 



TABLE V: Predictions for a few neutron-star observables. The various quantities are as follows: p t is the transition density 
from non-uniform to uniform neutron-rich matter, 7? is the radius of a 1.4 solar-mass neutron star, M max is the maximum mass, 
Purca i s the threshold density for the direct Urea process, and Murca is the minimum mass neutron star that may cool down 
by the direct Urea process. 



V. CONCLUSIONS 



The structure of neutron stars depends sensitively on the equation of state of neutron-rich matter in beta equilibrium. 
As such, nuclear effective interactions play a critical role in astrophysical applications. Given the large extrapolations 
in isospin and density that are required to compute the EOS, it is imperative to constrain the effective interaction. 
Motivated by two recent analyses that have provided simultaneous mass and radius information on neutron stars |16l 
117] . we have introduced a new effective interaction constrained by the ground state properties of finite nuclei, their 
collective excitations, the properties of both dilute and dense matter, and the structure of neutron stars. Thus, 
experimental, theoretical, and observational constraints have all been incorporated into the calibration procedure. 

Following the conclusions of Ref. [T7] — that suggest that the accurately calibrated FSU effective interaction is 
slightly stiff at intermediate densities but soft at high densities — we have used two empirical parameters to correct 
these shortcomings. The new effective interaction — dubbed "IU-FSU" — softens the EOS at intermediate densities by 
reducing the neutron skin thickness of 208 Pb and stiffens the EOS at high density by increasing the maximum neutron 
star mass relative to the FSU predictions. 

Existing relativistic mean-field interactions are flexible enough to reproduce, not only the charge density and binding 
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energy of closed shell nuclei, but also the neutron skin thickness of 208 Pb and the maximum neutron star mass (101141) . 
In particular, the density dependence of the symmetry energy is highly sensitive to the empirical parameter A v . Thus, 
an increase in the value of A v softens the EOS at intermediate densities and generates a neutron skin thickness in 
208 Pb of only R n — R p = 0.16 fm. For comparison, the NL3 and FSU effective interactions predict significantly larger 
neutron skins: 0.28 fm and 0.21 fm, respectively. To stiffen the EOS at high densities we rely on the empirical 
parameter £ [221 113 S3] ■ Including £ may be used to change the maximum neutron star mass by almost one solar 
mass without adversely affecting the bulk properties of nuclear matter around saturation density. Relative to the 
FSU model, the new IU-FSU interaction increases the maximum neutron star mass from 1.72 M to 1.94 Mq. 

As more accurate data becomes available, the coupling constants of the theory may need to be re-adjusted. The 
calibration procedure outlined here can accomplish this task rapidly and efficiently Ultimately, however, one would 
rely on an accurately calibration procedure to determine the best-fit parameters of the theory and to obtain correla- 
tions among the various observables. Nevertheless, as it stands now the new IU-FSU relativistic effective interaction 
reproduces: (a) the binding energies and charge radii of closed-shell nuclei, (b) various nuclear giant (monopole and 
dipolc) resonances, (c) the low-density behavior of pure neutron matter, (d) the high-density behavior of symmetric 
nuclear matter, and (e) the mass-radius relation of neutron stars. In the future, we plan to use the IU-FSU inter- 
action to generate a finite-temperature/finite-density equation of state for use in supernova and neutron-star merger 
simulations [42l 143] . 

This work was supported in part by DOE grants DE-FG02-87ER40365 and DE-FG05-92ER40750. 
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